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SUMMARY 


Two control approaches for a Large Gap Magnetic Suspension System (LGMSS) are investigated and 
numerical results are presented. The approaches are based on Linear Quadratic Regulator (LQR) control 
theory and Include a nonzero set point regulator with constant disturbance input and an integral feedback 
regulator. The LGMSS provides five degree-of-freedom control of a cylindrical suspended element which 
is composed of permanent magnet material. The magnetic actuators are air core electromagnets mounted 
in a planar array. 


INTRODUCTION 

This paper describes two control approaches for a Large Gap Magnetic Suspension System (LGMSS). 
The approaches are based on LQR control theory and include a nonzero set point regulator with constant 
disturbance input and an integral feedback regulator. The LGMSS is a conceptual design for a ground- 
based experiment which could be used to investigate the technology issues associated with magnetic 
suspension at large gaps, accurate suspended element control at large gaps, and accurate position sensing 
at large gaps. This technology would be applicable to future efforts which could range from magnetic 
suspension of wind tunnel models to advanced spacecraft experiment isolation and pointing systems. The 
LGMSS provides five degree-of-freedom control. The suspended element is a cylinder which is composed 
of permanent magnet material and the magnetic actuators are air core electromagnets mounted in a 
planar array. This co nfig uration is described and an analytical model developed in reference 1. Values of 
magnetic fields generated by the electromagnets at the location of the suspended core are presented and 
the method of calculation discussed in Appendix A which was written by Colin P . Britcher , Department 
of Mechanical Engineering and Mechanics, Old Dominion University, Norfolk, Virginia. All numerical 
results which are presented were obtained by using MATRIX x with SYSTEM -BUILD. 1 The MATRIXx 
implementation of the system is presented in Appendix B. Numerical results were obtained for positioning 
accuracy, with fixed gains, over large angular displacements about the yaw axis. 


LGMSS ANALYTICAL MODEL 

T his section presents an analytical model of the LGMSS. For a detailed development of the model see 
reference 1. The suspended element, or core, is a cylinder which is composed of permanent magnet 
material. The core is suspended over a planar array of five electromagnets mounted in a circular 
configuration. Figure 1 is a schematic representation of this system which shows the coordinate systems 
and initial alignment. A set of orthogonal x, y, z, body fixed axes defines the motion of the core with 
respect to inertial space. The core coordinate system is initially aligned with an orthogonal x, y, z system 
fixed in inertial space. A set of orthogonal X [ ,, y^, z& axes, also fixed in inertial space, define the location 
of the electromagnet array with respect to the x, y, z system. The x& and y& axes are parallel to the x and 
y axes respectively and the z f, and z axes are aligned. The centers of the two axis-systems are separated 
by the distance h . The angular acceleration of the core, in core coordinates, can be written as (see ref. 1) 

{ft} = (l// c ) (Vo!([W][T m ]{B}) + {T d }) (1) 

where I c is the core moment of inertia about the axes of symmetry (y and z), Vol is the volume of the 
core, M is the magnetization of the core, [M] is the skew symmetric cross product matrix, [‘ T m ] is the 
vector transformation matrix from inertial to core coordinates, {£} is the flux density produced by the 


^Use of names of products in this report does not constitute an official endorsement of such products, either expressed or 
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electromagnets, and {T d } represents external disturbance torques. A bar over a variable indicates that it 
is referenced to core coordinates. The translational acceleration of the core, in core coordinates, can be 
written as 

= 0 /»«<•) (vol ([r„,][dB][r,„]- 1 { J V/}) + ( 2 ) 

where m c is the mass of the core, \0B] is a matrix of the gradients of B , and { F d } represents external 
disturbance forces in core coordinates. Under the assumptions of reference 1, {/?} can be written as 

{5} = (l//max) [A/j]{/} (3) 

where / max is the maximum coil current, [A'g] is a 3 x 5 matrix whose elements represent the values of 
{B} produced by a corresponding coil driven by the maximum current and the coil currents are 

{I} T = [h I 2 h h h\ (4) 

The gradients can be put in the same form by arranging the elements of [OB] as a column vector. This 
results in 

{dB} = (lfI max )[K im }{I} (5) 

where {cU?} is a nine element column vector containing the gradients of {£} , and [Kob] is a 9 x 5 matrix 
whose elements represent the values of {OB} produced by a corresponding coil driven by the maximum 
current. The simplified notation from reference 1, where B tj represents OBjdj , will be used for gradients. 
Each element of { dB }, for example B xx , can be written in the form 

B xx = (1/Anax) (6) 

where [A XX J is a 1 x 5 matrix containing values of B xx produced by a corresponding coil driven by the 
maximum current. A block diagram of the system is shown in figure 2. This model is nonlinear because 
of the coordinate transformations and is of the form 

x = f{x, u) (7) 

where x is given by 

x T — [ft,, ft r V x Vj, K 0 y 0 : x y z\ (8) 

and the input u is given by 

u 1 - [/i I 2 h h h J (9) 

The states and $ z in (8) are the pitch and yaw angles of the core respectively. Analyses and 
simulations of the model were performed using MATRIX y with SYSTEM-BUILD. A description of the 
MATRIXx /SYSTEM -BUILD implementation is presented in Appendix B. LGMSS model parameters are 
presented in Table Al. 


CONTROLLER DESIGN APPROACH 

The control approach which was investigated was to linearize the system about a nominal operating 
point x a , u 0 and then use the resulting linear model to calculate feedback gains to stabilize and control 
the core about this point. The equations for computing u v for a given x 0 are presented in the appendix 
of reference 1. The feedback gains for the linearized system were calculated using LQR theory. In this 
approach the system is put in state variable form 

x = Ax + Bu (10) 


and a set of full state feedback gains arc calculated that minimizes the quadratic cost function 



where the weighting matrices, Q and /?, are symmetric and positive-scmidcfinitc and positive-definite, 
respectively. The resulting control law is 

a — —KjfX ( 12 ) 

where 

K, { = R~ ] B r P (13) 

and P is a symmetric positive-definite matrix satisfying the matrix Riccati equation (see ref. 2) 

PA 4- A r P + Q- PBlt - l B r P - 0 (14) 

The operation of the LGMSS requires that the control system provide variable set point control and 
also that it compensate for a constant disturbance (the weight of the core). Two feedback controller 
configurations which meet these requirements were investigated and arc described below. 

Nonzero set point regulator with constant disturbance.- With a constant disturbance, V 0 , the system 


under consideration becomes 

x — Ax -1- Bu + V„ 

(15) 

The controlled variables arc assumed to be 

2 = Dx 

(16) 

For the LGMSS, z is defined as 

z 7 = [By 6 : x y z\ 

(17) 

For a set point of z Q , a constant input u Q 

must be found to hold the states at x Q so that 



z 0 — Dxq 

(18) 


From (15) above, x Q and u a must be related by 


0 = -}- Bu 0 -f- V Q (19) 

Next, define the shifted input, shifted state, and shifted controlled variable, respectively, as 

u = u - ii 0 (20) 

x — x - x 0 (21) 

z' = 2 - (22) 

Solving equations (20)~(22) for u, x , and z, substituting the results into (15) and (16) and using (18) and 
(19), it can be shown that the shifted variables satisfy the equations 

x f = Ax' + Bu' + Vo (23) 

and 

z' - Dx’ (24) 

Assume the shifted regulator problem has a steady-state solution in the form of the control law 

u! = —Fx (25) 


that minimizes the quadratic cost function 


where F is a matrix of constant feedback gains. Prom (25) and (20) -(22) 


u - -Fx + u 0 + Fx 0 (27) 

which has the form 

u = -Fx + u' a (28) 

where the constant vector u' 0 is to be determined so that in steady-state the controlled variable z assumes 
the given value z 0 . Substituting (28) into (15) gives 

x = (A — BF)x + Bu 0 + V 0 (29) 

Since the closed-loop system is assumed stable, as t — > oo the states approach x 0 which satisfies 

0 = Ax 0 -f Bu 0 + V 0 (30) 

where 

A = A - BF (31) 

For a stable system A is nonsingular and x Q becomes 

, = {-A)- l Bu’ 0 + (-A)- 1 V 0 (32) 

z 0 then becomes 

*o = Di-Ay'Bu'e + Dl-Ay'Vo (33) 

and solving for u' 0 results in 

u' 0 = [D(-A)~ l B)- l z 0 - [Di-Ar^^Di-AyWa (34) 

A block diagram of this controller is shown in figure 3. 

Integral feedback regulator .- Another controller configuration which was investigated is one which uses 
integral feedback. Integral feedback is often used in systems subjected to constant disturbances in order 
to minimize steady-state errors (ref. 2). To develop the integral feedback controller equations, assume the 
system of (15) with the controlled variables defined by (16). To the system of (15) add the integral states 
q , defined by 

q = z (35) 

where z' is defined by (22). A performance function can now be defined as 


roc 

J = / ( z t Qz + (p Qq 4- u T Ru)dt 

JO 

Assume that the time-invariant control law 


(36) 


u = -Fix - F 2 q (37) 

can be found which stabilizes the augmented system described by (15), (16), and (35), and which minimizes 
(36). Since the presence of the constant disturbance vector in (15) has no effect on the asymptotic stability 
of the system 

lim q = 0 (38) 

t — ♦ oo 

and from (35) 
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lim z = O 
t—>oo 


(39) 



This means that the control system with the control law (37) has the property that the error in the 
controlled variables, due to a constant input disturbance, eventually goes to zero. A block diagram of the 
integral feedback regulator is shown in figure 4. 

Initial conditions .- The initial conditions on current for the nonzero set point regulator with a given 
x 0 and z 0 can be determined from equation (34). The initial conditions on the integrator outputs, for 
the integral feedback regulator, can be determined from (15) and (37). In equilibrium, assuming a stable 
system, (15) reduces to (19). Substituting (37) into (19) results in 


0 = Ax 0 - BF 2 q 0 + V 0 

(40) 

where A = A — BF\. Solving for x 0 results in 


x 0 = A~ l BF 2 q 0 - A~ l V 0 

(41) 

and z 0 becomes 

z 0 = DA- 1 BF 2 q 0 - DA-W 0 

(42) 

Finally, q 0 becomes 

q a = ( DA~ l BF 2 ) _1 Zo + ( DA- l BF 2 y l DA~ l V 0 

(43) 


Discrete model .- The equations presented above are for a continuous system and assume a perfect 
sensor with zero delay. In order to investigate the effects of a delay introduced by a finite sensor processing 
time, a worst case sampling time of 20 samples per second was assumed. The continuous system equations 
were discretized and discrete system gains were calculated using this sample time. The sensor delay time 
was modeled as a zero order hold (ZOH) on the output of the nonlinear open-loop system for simulation 
purposes. A block diagram showing the placement of the ZOH for the integral feedback regulator is shown 
in figure 5. 


ANALYSIS AND SIMULATION RESULTS AND DISCUSSION 

This section presents the analysis and simulation results for the nonzero setpoint and integral feedback 
regulator which were obtained by using the MATRIXx implementation described in Appendix B. The 
results are for positioning accuracy of the the two approaches, with fixed gains, over large angular 
displacements about the yaw axis. The approach was to define, first of all, the range over which a 
set of fixed gains was stable for each regulator type. This was accomplished by obtaining linearized 
models of the open-loop system at intervals of two degrees yaw. Feedback gains were computed for a 
given reference angle and these gains were then used with the linearized models to obtain closed-loop 
eigenvalues over a given range of yaw angles approaching the limits of stability. The full nonlinear model 
was then simulated using feedback gains calculated at a given yaw angle. Inputs to the system consisted of 
position commands in one degree increments. For more detail on the design model used, see Appendix B. 


Continuous System 

Nonzero set point regulator.- A plot of the eigenvalues of the nonzero set point regulator with feedback 
gains calculated for the model linearized around zero yaw angle, as the core is rotated through increments 
of two degrees, is shown in figure 6. The low frequency eigenvalues are shown on an expanded scale 
in figure 7. The point at which the eigenvalues cross over the imaginary axis is at a yaw angle of 
approximately 42 degrees. A plot of the eigenvalues, with the same feedback gains, as the model is 
rotated from zero degrees to —42 degrees is shown in figure 8. The low frequency eigenvalues are shown 
on an expanded scale in figure 9. As expected, because of symmetry, these eigenvalue plots are identical 
with the previous plots. Further calculations also verified that this pattern repeats every 36 degrees. That 
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is, fixed gains calculated at zero, 36, 72, ect. are stable over a range of plus and minus 42 degrees around 
the angle for which the gains are calculated. The system performance was evaluated by commanding yaw 
position over a range of 18 degrees in increments of one degree. Time between increments was 20 seconds. 
Figure 10 is a plot of the input command and figure 11 is a plot of the actual yaw angle. Figure 12 is a 
plot of the difference between command and actual angle. Figures 13 through 16 show the response of the 
pitch, x, y, and z axes. As can be seen, the errors build up as the yaw angle increases. This is the main 
disadvantage of the nonzero set point controller for this particular application. The input which minimizes 
the position error for a given command and constant input disturbance is given by equation (34) and is a 
function of A which is in turn a function of A. As evident from the block diagram of the analytical model 
in figure 2, A is a function of core position. The system response to a command over a negative 18-degree 
range was essentially the same and is not shown. The system response about other nominal gain points, 
at 36-degree intervals in yaw, was also essentially the same. 

Integral feedback regulator .- A plot of the eigenvalues for the integral feedback regulator with fixed 
gains calculated for the model linearized around zero yaw angle is shown in figure 17. The low frequency 
eigenvalues are shown on an expanded scale in figure 18. The symmetry properties for the nonzero set 
point regulator discussed above were found to apply to the integral feedback regulator also. The response 
of the system with the integral feedback controller with fixed gains calculated at zero degrees, to the 
co mm and input shown in figure 19, is shown in figure 20. The yaw error is shown in figure 21 and the 
response of the pitch, x, y, and z axes are shown in figures 22 through 25. As can be seen from these 
figures, and as predicted by equation (39), after an initial transient the errors all go to zero. As was 
the case with the nonzero set point regulator discussed above, system response about other nominal gain 
points (at 36-degree intervals in yaw) was essentially the same. 


Discrete System 

Integral feedback regulator .- After evaluating the performance of the nonzero set point regulator with 
a continuous system model, a decision was made to drop this approach from further consideration and 
concentrate on the integral feedback regulator for the discrete system evaluation. A plot of the eigenvalues 
of the discrete integral feedback regulator with feedback gains calculated for the model lin earized around 
zero yaw angle, as the core is rotated through increments of two degrees, is shown in figure 26. The 
eigenvalues shown have been mapped from the z-plane into the s-plane. The low frequency eigenvalues 
are shown on an expanded scale in figure 27. As can be seen from the figures, the range of stable operation 
is reduced from 42 degrees for the continuous system to approximately 36 degrees for the discrete system. 
The discrete system was symmetric about points 36 degrees apart in yaw angle as was the continuous 
system. The response of the discrete integral feedback controller with fixed gains, calculated for the model 
linearized around zero degrees, to the command input shown in figure 28 is shown in figure 29. The yaw 
error is shown in figure 30 and the response of the pitch, x, y, and z axes are shown in figures 31 through 
34. The transients in the response of the discrete system were of higher magnitude than the continuous 
system but not by a significant amount. Other than magnitudes, the responses were similar with the 
errors returning to zero after initial transients. 


CONCLUDING REMARKS 

Two control approaches for a Large Gap Magnetic Suspension System (LGMSS) have been investigated 
and numerical results presented. The control approaches investigated were a nonzero set point regulator 
with constant disturbance input and an integral feedback regulator. The LGMSS provides five degrees 
of freedom control of a suspended element which is a cylinder composed of permanent magnet material. 
The magnetic actuators are air core electromagnets mounted in a planar array. Results were obtained 
for both continuous system and discrete system models. The continuous system model assumed a perfect 
position sensor with no delays and the discrete system model assumed a worst case delay introduced 
by the sensor of fifty milliseconds. All analyses and simulations were performed using MATRIXx with 
System_BUILD. The numerical results are for positioning accuracy, with fixed gains, over large angular 



displacements about the yaw axis. These results indicate that for the LGMSS investigated, the integral 
feedback regulator provides the best performance. Also, for control of yaw over a range of 360 degrees, 
scheduled gains will be required. Based on the numerical results, control about a nominal operating point 
is possible over a range of plus or minus 18 degrees. Although the system is stable over a greater range, 
transient errors become larger as the yaw angle is increased and plus or minus 18 degrees appears to 
be a reasonable limit. Gains can be scheduled over a smaller interval based on particular requirements. 
Control of pitch over large angles was not addressed in this investigation and should be the subject of a 
future effort. 

The results presented in this paper, for the analytical model utilized, indicate that delays introduced by 
the position sensor, for the worst case value investigated, can be adequately compensated. Future efforts 
should include investigations of the effects on performance of: 1) different output feedback approaches, 
2) changes in {B} and [#B] with core displacement, 3) power supply and coil dynamics, and 4) sensor 
hardware characteristics. 
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APPENDIX A 


CALCULATION OF ELECTROMAGNET FIELDS AND GRADIENTS 

By 

Colin P. Britcher 

Department of Mechanical Engineering and Mechanics 
Old Dominion University 
Norfolk, Virginia 

Magnetic fields external to air-cored electromagnets of simple geometry can be calculated accurately 
by the use of line conductor models. A computer program FORCE (ref. 5) represents such electromagnets 
by an assembly of straight line conductor elements. Fields and field gradients produced by each element 
at any point, not coincident with the element, may then be calculated exactly using the Biot Savart law. 
Fields and field gradient components are then simply summed over all conductor elements. This approach 
is illustrated in figure Al. 

This procedure is accurate if the number of conductor elements is large enough to properly represent the 
three-dimensionality of the electromagnet yet small enough to avoid excessive accumulation of numerical 
(rounding) errors. Accuracy generally improves with distance from the electromagnet as the effects of 
the discontinuous representation of current are reduced. Recalculation with greater or lesser numbers of 
elements is a straightforward and effective check on the reliability of the solution. 


Calculations using FORCE were made for the electromagnet configuration used in the earlier stages 
of the design study reported in reference 6. The parameters for this system are summarized in table Al 
and a schematic representation is given in figure A2. Calculation results are presented in table A2. It 
should be noted that detail adjustments to the configuration and to the design suspension height were 
made later in the design study mentioned above but do not substantially affect the results shown here or 
the conclusions of this report. 


TABLE Al.- LGMSS MODEL PARAMETERS 


Core diameter, m 

Core length, m 

Core mass (m c ), kg 

Core inertia ( I c ), kg-m 2 

Core volume (Vol), m 3 

Core magnetization (A A/m . . . . 

Core suspension height (h), m 

Electromagnet outer radius, m 

Electromagnet inner radius, m 

Electromagnet height, m 

Maximum electromagnet current (/ m ax)> A 


. . 0.1016 

. . 0.3048 

. . 22.76 
... 0.6 
2.5(10)" 3 
9.5493(10) 5 
. 0.9144 

. . 0.386 
. . 0.173 
. . 0.493 
. . 559.5 
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TABLE A2. - ELECTROMAGNET FIELDS AND GRADIENTS 


Electromagnet 

Bx 

By 

By, 

B X x 

B\Y 

B xz 

Byy 

Byz 

B zz 

(Tesla) 

(Tesla/meter) 

1 


o 

■ 


0 

+.0544 




2 

E 

■ 


+ .0296 

-.0129 

+.0168 


+.0518 


3 



tSkm 

+.0051 

+.0209 

-.0440 


ami 


4 

+.0192 

E 

+.0218 

+.0051 

-.0209 

-.0440 




5 

-.0073 

+.0225 

+.0218 

+.0296 

+ .0129 

+.0168 

-.0059 

-.0518 

HSa 


Fields correspond to maximum rated current for each 


electromagnet (559.5 Amp.). 
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0.9144m 



Figure A1 .-Geometry used for magnetic field calculations. 



Figure A2.-Representation of electromagnets by line conductor elements. 


APPENDIX B 


MATRIX x IMPLEMENTATION 

qwrwMmm n lmp] * ment ™Z the LGMSS plant model and control approaches, using MATRIXx and 

, . -BUILD, involves the construction of block diagram representations of the system (figure Bl) 
and the controllers (figures B2, B3, and B4). Since version 6.0 of MATRIX x /SYSTEM JBUILD was used 
to construct the models, it was necessary to break the implementation into nested super blocks containing 
a maximum of 6 blocks each (for more detailed information on MATRIX x /SYSTEM _BUILD see refs. 3 
an 4). Figure B5 shows the actual super blocks of the implementation. Each of the blocks in figure Bl 
can e identified as a block within some level of the nested super blocks shown in figure B5 For the 
sake of brevity and clarity only the unnested versions of the SYSTEM .BUILD models will be shown and 
discussed from this point on. 

There is a correspondence between blocks in the block diagram of the system (figure 2) and blocks 
or groups of blocks, in the SYSTEM J3UILD implementation of figure Bl. For instance, the leftmost 
block m figure 2 is modeled by the super block labeled FIELDS in figure Bl. The expression (1/J max ) is 
represented by the gain block labeled I MAX in which each of the 5 outputs is calculated by multiplying the 
corresponding input by a constant gain factor of (l// max ). The K matrix multiplication is implemented 
by the State Space block labeled XY Z which multiplies the vector of 5 currents by a 5 x 12 matrix to 
calculate the vector of 12 outputs corresponding to the 3 field components (B x , B y , B z ) and the 9 gradients 

( Bxx ,_c)r?r?r?DDr>r>\ a ° 


xy 


> Bxz > Byxi Byy, By Z , B ZX , B Z y 




One block found in figures Bl and B5 has no corresponding element in the system block diagram of 
figure 2. This block, labeled A GRAV outputs a vector [0 0 g] which is passed through a transformation 
from inertial to body coordinates (block TM) and summed with the linear accelerations to model the 
effects of gravity on the core. 


In order to set up for calculation of the LQR gains, the SYSTEM .BUILD ANALYZE option is invoked 
on the super block OPEN (fig. B6). This super block contains a super block labeled PLANT, which is 
the level super block in the plant model of figure B5 (and represents the whole of fig. Bl). Super 

block OPEN also has a step input block labeled I OFF which sets the initial currents about which the 
nonlinear system will be linearized and a gain block labeled CONTIN. The CONTIN block has unity gain 
and acts as a place holder for a zero order hold super block used in discrete versions of the system as will 
be discussed below. The yaw angle ( 9 Z ) for linearization is established by setting initial conditions on the 
corresponding integrator prior to analysis of the model. 


TT J he . A ^ ALYZE operation returns from the SYSTEM .BUILD module to MATRIX x where the 
earize function is invoked to calculate the linearized state space representation of the nonlinear system 
represented by PLANT. LIN returns an S matrix and NS, the number of states in S (10 for this model). 
I he state space matrix S can be partitioned into 


- 

A 

‘ 

- C 

_ 

D 


(Bl) 


where A and B are the matrices in the state variable representation of the system in equation 10. The 
matrix A will be of dimension NS x NS, B will be dimensioned NS x NI, C will be NO x NS and 
^ wi |! ke NO x NI; where NI is the number of inputs (5 currents) and NO is the number of outputs 
((A, I , Z) positions and velocities, Y and Z rotation angles and rotational velocities). The MATRIXx 
function SPLIT (5, N S) accomplishes this partitioning and can return A and B only or all four matrices 
For the problem being considered only A and B are of interest since C turns out to be an identity matrix 



and D is a zero matrix. The output equation 


y = Cx + Du (B2) 

reduces to 

V = x (B3) 

since the outputs are the position and velocity variables corresponding to the states, and outputs are not 
directly related to inputs. For the model linearized around zero, the A matrix becomes 



162.17 

0 

0 

0 

0 

5x5 i 

0 

162.17 

0 

0 

0 

[0) 

-9.8066 

0 

0 

0 

0 


0 

0 

0 

0 

0 


-0.0147 

-0.0049 

0 

0 

0 

5x5 

. r/i i 


5x5 

[01 





and the B matrix becomes 


' -0.155 

-0.155 

-0.155 

-0.155 

-0.155' 

0 

-0.16 

-0.098 

0.098 

0.16 

-0.0019 

0.0055 

0.001 

0.001 

0.0055 

0 

-0.0024 

0.0039 

-0.0039 

0.0024 

0.102 

0.0031 

-0.0082 

0.0082 

0.0031 



5x5 



. 


[01 




This model has four nonzero eigenvalues: -12.7346, -12.7346, 12.7346, and 12.7346. 


(B4) 


(B5) 


A SYSTEM J3UILD implementation of the linearized model was constructed for use in testing and is 
shown in figure B7. The full nonlinear model was employed for all simulations shown in this document. 


The nonzero set point regulator gains can then be calculated using the linearized model represented 
by A and B, and the Q and R matrices from the quadratic cost function of equation 26. The Q and R 
matrices used for gain calculation were 

Q = [l(10) 6 J (B6) 

and 

(B7) 

where Q is a 10 x 10 diagonal matrix and R is a 5 x 5 identity matrix. The MATRIXx function 
REGULATOR calculates the LQR gain matrix F and the corresponding eigenvalues, EV. As a baseline, 
the closed-loop eigenvalues for the system with the model linearized around zero and gains calculated 
using (B6) and (B7) are presented in figures B8 and B9. The nonzero set point matrix HCI can then be 
calculated. HCI is given by 

HCI = [£>(-l4) _1 .B] _1 (B8) 

where D is defined in equation 16. For this implementation the only nonzero elements in D are those 
corresponding to A r , Y, Z, 9 y , and 8 : . The value of each nonzero element is one. The SYSTEM _BUILD 
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implementation of the nonzero set point controller is shown in figure B2 and corresponds to the 
conventional block diagram of figure 3. Currents as well as position and rotation variables are brought to 

the output. Inputs and outputs are scaled to appear in inches and degrees while internal quantities are 
meters and radians. 

The integral feedback continuous LQR gains were calculated similarly by augmenting the linear state 
space model, represented by .4 and B , with 5 new states corresponding to the feedback integrators: 



10x10 

5x5 

A' = 

[-'1] 

[0] 


5x5 

5x5 

5x5 


10] 

r/j 

[0] . 


rlOxfi n 



5x5 


L [0] J 


(B9) 


(BIO) 


The MATRIXx function REGULATOR was again used to calculate LQR gains. The closed-loop 
eigenvalues for the integral feedback system with the model linearized around zero and gains calculated 
using (B6), extended to a 15 x 15 diagonal matrix, and (B7) are presented in figures BIO and Bll. 
Figure B3 shows the SYSTEM-BUILD model of the continuous integral feedback controller of figure 4. 
The delay due to the sensor was modeled by using the MATRIX x DISCRETIZE function which takes as 
parameters a state space matrix, 5, its order, NS, and a delay TAU which was specified as 0.05 seconds, 
and returns a discrete state space matrix SD. A state space matrix for the augmented system, S' 
was constructed using A' and B' above as well as dummy C' and D' . The discretized version of’this 
matrix, SD, is then SPLIT and the resulting AD and BD system matrices used with MATRIX x function 
DREGULATOR, similar to REGULATOR discussed above, to calculated LQR feedback gains for the 
discrete version of the system. The closed-loop eigenvalues for the discrete system with the model 
linearized around zero and gains calculated using (B6), extended to a 15 x 15 diagonal matrix, and 
(B7) are presented in figures B12 and B13. The eigenvalues shown have been mapped from the 2 -plane 
into the s-plane. Figure B14 shows the SYSTEM _BUILD model of the discrete integral feedback controller 
of figure 5. The discrete super block ZOH is the sensor model. It replaces the CONTIN block and acts 
as a Zero Order Hold. ZOH is the only discrete super block used in this otherwise continuous model. 
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Figure B1 .-Nonlinear open-loop system. 






































Figure B3. -Integral feedback regulator 

















































Figure B5 .-Nested super block model of nonlinear open-loop system 














Figure B7. -Linearized open-loop system. 
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Figure BIO. -Closed-loop eigenvalues for the integral feedback 
system with the model linearized around zero. 







- 



X 

X- 





- 












Figure B1 4. -Super Block implementation of intergral feedback regulator 
with zero order hold. 
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Figure 1 .-Initial coordinate system alignment for large gap magnetic 
suspension system. 
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Figure 2.-Block diagram of analytical model of large gap magnetic 
suspension system. 
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Figure 3.-Nonzero set point regulator. 
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Figure lO.-Input to nonzero set point regulator with feedback gains 
calculated at zero yaw angle. 
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Figure 1 1 .-Yaw angle response to command input of figure 1 0. 
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Figure 1 2.-Difference between command input of figure 1 0 and 
response of figure 11. 
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Figure 1 3.-Pitch angle response to command input of figure 1 0. 
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Figure 15.-Y axis translational response to command input of 
figure 10. 
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Figure 16.-Z axis translational response to command input of figure 10. 
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Figure 18.-Low frequency eigenvalues of figure 17. 
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Figure 21 .-Difference between command input of figure 19 and response of 
figure 20. 
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Figure 22.-Pitch angle response to command input of figure 19 
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Figure 23.-X axis translational response to command input of figure 19. 
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Figure 24.-Y axis translational response to command input of figure 1 9. 
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Figure 27,-Low frequency eigenvalues of figure 
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Figure 31 .-Pitch angle response to command input of figure 28. 
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Figure 33 .-Y axis translational response to command input of figure 28. 





IWV5A 

Naionjy Amman and Report Documentation Page 

SpacaAdnwMfriKr 


1 Report No. 2. Government Accession No 

NASA TM- 101 606 


•1. Title and Subtitle 

An LQR Controller Design Approach for a Large Gap Magnetic Suspension 
System (LGMSS) 


3 Reaptenrs Catalog No. 


5. Report Date 
July 1990 


6. Performing Organization Code 


t Author(s) 

Nelson J. Groom and Philip R. Schaffner 


8. Performing Organization Report No. 


9 Performing Organization Name and Address 
NASA Langley Research Center 
Hampton, VA 23665-5225 


10. Work Unit No. 
505-66-91-02 


11. Contract or Grant No. 


12 Sponsonng Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


13 Type of Report and Period Covered 

Technical Memorandum 

14. Sponsonng Agency Code 


15. Supplementary Notes 


16. Abstract 

Two control approaches for a Large Gap Magnetic Suspension System (LGMSS) are investigated and numerical 
results are presented. The approaches are based on Linear Quadratic Regulator (LQR) control theory and include a 
nonzero setpoint regulator with constant disturbance input and an integral feedback regulator. The LGMSS which is 
considered provides five degree of freedom control. The suspended element is a cylinder which is composed of 
permanent magnet material and the magnetic actuators are air core electromagnets mounted in a planar array. 


17 Key Words (Suggested by Author(s)) 

Magnetic Suspension; Magnetic Levitation; 
Pointing System; Magnetic Suspension Model; 
Magnetic Suspension Control System 


18. Distribution Statement 


Unclassified - Unlimited 
Subject Category 31 


19 Security Classtf. (of this report) 

Unclassified 


20. Security Classif. (of this page) 

Unclassified 


21. No. of pages 

64 


22. Price 


NASA FORM 1626 OCT B6 For sale by the National Technical Information Service, Springfield, VA 22161-2171 










